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Abstract. - We show that a broad class of quantum critical points can be stable against locally 
correlated disorder even if they are unstable against uncorrelated disorder. Although this result 
seemingly contradicts the Harris criterion, it follows naturally from the absence of a random- 
mass term in the associated order-parameter field theory. We illustrate the general concept with 
explicit calculations for quantum spin-chain models. Instead of the infinite-randomness physics 
induced by uncorrelated disorder, we find that weak locally correlated disorder is irrelevant. For 
larger disorder, we find a line of critical points with unusual properties such as an increase of the 
entanglement entropy with the disorder strength. We also propose experimental realizations in 
the context of quantum magnetism and cold-atom physics. 



Introduction — The effects of quenched disorder in con- 
L^jdcnscd matter are various and constitute an interesting 
and important field of research. For instance, randomness 
can change the universality class of a critical point and 
even originate novel phases. A paradigmatic effect of dis- 
qq order near a phase transition is known as "random mass:" 
H As there is no translational invariance, different regions of 
C>the same system can be at different distances from criti- 
^— H cality 1 . Thus, disorder fluctuations induce large and rare 
regions which can be "locally" in one phase while the bulk 
( ^ is in the other one. This yields the so-called Griffiths sin- 
. . gularities [1-3]. 

^ • Using the random-mass concept, Harris [4] formulated 
a simple criterion for the relevance of disorder at continu- 
_^ ous phase transitions: If dv < 2 (with d being the spatial 
C3 dimension and v the clean correlation- length exponent), 
then the clean critical behavior is destabilized by weak 
disorder. This is the famous Harris criterion. It applies 
to the case of uncorrelated disorder. Spatial correlations 
among the random masses modify the criterion. When 
such correlations decay as x~ a with distance x, the Harris 
criterion reads min{d, a]v < 2 [5]. This implies that corre- 
lations are relevant only when they decay slower than x~ d . 
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*In field-theory language, it is as if the mass of the field were 
random in space. For thermal phase transitions, this idea is also 
known as random- T c . 



Interestingly, when the disorder correlations are relevant, 
the Griffiths singularities are also enhanced [6]. From the 
above arguments, one may expect that local (i.e., short- 
range) disorder correlations do not change the relevance 
or irrelevance of the disorder. 

In this Letter, we show that this is not always true. 
Specifically, we demonstrate that certain types of local 
correlations can render the disorder perturbatively irrel- 
evant even though dv < 2, thereby stabilizing the clean 
critical point against weak disorder. Although this result 
appears to contradict the Harris criterion, it arises natu- 
rally when disorder correlations make the mass term spa- 
tially uniform. If the strength of the correlated disorder 
is increased beyond a critical value, where perturbative 
methods cannot be used anymore, a line of finite-disorder 
critical points (tuned by disorder strength) appears. 
Field theory — Before turning to a specific microscopic 
model, let us consider the relevance of disorder within a 
general field theory framework. Consider a (clean) Eu- 
clidean quantum- field theory given by the action So[</>] = 
J drdyiCo [(f)] (x represents position in d dimensions and 
r is imaginary time). The action Sq is perturbed by a 
disorder term 5dis = J2o I drdx Ae>(x)0(x, t). Here, the 
sum is over all possible operators O to which disorder can 
couple, and Xq (x) is a random variable of zero mean. (As 
we are considering quenched disorder, \q depends only 
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on the spatial coordinates.) Using the replica trick and 
performing a tree-level renormalization-group (RG) cal- 
culation (see, e.g., Ref. [7]), the relevance of Ao(x) can 
be gauged by computing the scale dimension [wq] of the 
second moment wq of Ae>(x) at the clean fixed point. If 
[wq] = d+2zo — 2[0] > 0, where [O] is the scale dimension 
of O and z$ is the dynamical exponent of the clean theory 
So, then disorder is pcrturbatively relevant. This can be 
viewed as a generalization of the Harris criterion. 

Applying these ideas to a 4 order-parameter theory, 
£o = r\4>\ 2 + u\(f)\ 4 + \d x (j)\ 2 + \d T cj)\ 2 , we find that quenched 
disorder can be relevant at tree level only if it couples 
to the mass term \<fi\ 2 . 2 (We exclude random fields that 
locally break the symmetry.) This follows from the fact 
that [tuuia] = —d + 2/v, which recovers the original Har- 
ris criterion, while disorder coupling to higher order and 
gradient terms leads to negative scale dimensions and is 
thus irrelevant. Consequently, if the macroscopic random 
variables do not produce a random mass term, disorder 
effects are strongly suppressed. 

Transverse field Ising chain — A tantalizing example in 
which the above scenario actually completely removes the 
random mass and thus stabilizes the clean critical point 
even though the inequality dv < 2 is fulfilled is the ferro- 
magnetic quantum phase transition of the ID transverse- 
field Ising model. Its Hamiltonian reads 
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where Jj's are the interactions, hi's are the transverse 
fields, and af and of are Pauli matrices. The correlation- 
length exponent of the clean (J, = J, hi = h) critical point 
is v = 1, and the Harris criterion predicts that weak dis- 
order in the hi and Ji is relevant. In agreement with this, 
uncorrelated disorder has dramatic effects [2,8,9]. The 
clean critical point is unstable against weak disorder and 
the RG flows towards an exotic infinite-randomness criti- 
cal point (IRCP) whose dynamics is so slow that the dy- 
namical exponent z is formally infinite. Surrounding the 
transition, there are gapless quantum Griffiths phases in 
which the dynamical exponent z can be arbitrarily large 
and the correlation length £ is finite. Moreover, the aver- 
age entanglement entropy S [10, 11] of a subsystem with 
length I embedded in the bulk diverges at the IRCP as 
S{£) ~ (c e ff/3)ln£, where c c s = (1/2) ln2 is often called 
the effective central charge [12,13]. The properties of this 
model are summarized in Fig. 1(b). 

The order-parameter field theory of this transition can 
be obtained from the Hamiltonian (1) via a Trotter-Suzuki 
decomposition of the partition function and a Hubbard- 
Stratanovich transformation of the resulting path integral. 
After taking the continuum limit, the resulting (1 + 1)- 
dimcnsional action takes the form of a disordered cj> 4 the- 
ory (see above). The (local) coefficient of the mass term 



2 If higher-order terms are disordered, a random mass may be cre- 
ated at loop level, but for weak disorder it will be strongly suppressed 
compared to the generic case. 
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Fig. 1: (Color online) (a) Disorder averaged spectral gap A vs. 
distance from criticality S for different disorder parameters ho 
with and without disorder correlations (sizes up to 2 15 sites, 
averaged over 1 000 disorder realizations), (b) RG flow diagram 
for uncorrelated disorder, the disorder strength renormalizes to 
infinity, (c) Flow diagram for correlated disorder (hi = e s Ji), 
the clean fixed point is stable. The line of fixed points for 
strong disorder can be only accessed for ho = (see text). 



\4>\ 2 must be a function of the ratios hi/ Ji because these 
are the only two energy scales in the problem. In fact, it 
has been shown [14] that the critical point of (1) occurs ex- 
actly when Yl Ji — II hi, suggesting that Si — ]n(hi/Ji) is 
an appropriate measure of the local distance from critical- 
ity 3 . Consequently, it is clear that a simple local correla- 
tion between hi and Ji, say, hi = e s Ji, is sufficient to make 
the mass term uniform rather than random. The disorder 
of the higher order and gradient terms of the renormalized 
action cannot be inferred from this argument, generically 
they remain random because they do not depend on the 
combinations hi/ Ji only 4 . Although translational sym- 
metry is broken, the system is locally at the same distance 
from criticality everywhere. 

Numerical results — To check the prediction that such cor- 
related disorder is irrelevant at the clean critical point, 
we mapped the Hamiltonian (1) onto free fermions us- 
ing the Jordan- Wigncr transformation [15,16]. We then 
performed an exact-diagonalization study contrasting the 
cases of uncorrelated and correlated disorder. The absence 
of Griffiths phases can be verified by analyzing the gap A 
in the excitation spectrum as a function of the distance 
from criticality 5 = [£ f ] av = [ln(^/ Jj)] av [see Fig. 1(a)]. 
Here, the fields are drawn from a uniform (box) distribu- 
tion in which ho < h < 1. Hence, ho parameterizes the dis- 
order strength. For correlated disorder, we set hi = e* 5 Ji, 
while for uncorrelated disorder, the couplings arc indepen- 
dent random variables drawn from a uniform distribution 



3 At first glance, this definition seems to break the symmetry be- 
tween left and right neighbors in (1). However, the concept of a 
distance from criticality is well defined only for regions large com- 
pared to the lattice constant. After averaging <5; over such a region, 
the discrepancy vanishes. 

4 It is important to distinguish the bare and renormalized actions. 
Pfeuty's criticality condition [16] guarantees that the renormalized 
theory does not have a random mass term. This does not imply that 
this term must vanish in any bare theory. 
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such that, for the sake of comparison, h^e^ 5 < Ji < e~ 5 . 
Figure 1(a) clearly shows that the gap vanishes only at 
criticality (S = 0) for correlated disorder, signalling the 
absence of quantum Griffiths phases. In contrast, the gap 
closes before criticality is reached when the disorder is un- 
corrected, a hallmark of quantum Griffiths singularities. 

The above uniform distributions of fields and couplings 
represent weak or moderate disorder (as long as kg ^ 0). 
To study the fate of the critical point for stronger cor- 
related disorder (and, therefore, study nonperturbative 
effects of disorder not accomplished by the conventional 
ficld-thcorctic analysis), wc now consider a family of gap- 
less power-law distributions 
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V{h) = (l/D) h- 1+1/D , with < h< 1, 



(2) 



with the bonds given by Ji = e~ s hj as before. D param- 
eterizes the disorder strength. In this case, the system is 
gapless for any 5; the absence of a gap simply follows from 
the low-energy tail of the distribution V . 

We first focus on the correlation length £ as a function 
of S which is computed by fitting the spin-spin correlation 



function 
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Figure 2(a) shows that 



£ takes its clean value regardless of disorder strength D. 
The correlation length exponent is thus v — I'dean = 1 
for all D, i.e., it violates the inequality dv > 2 [17]. This 
implies that the spin correlations are governed by the clean 
physics analogous to gapped quantum magnets doped with 
non- magnetic impurities [18]. We now explore the critical 
point Ji = hi, upon increasing the disorder strength D. 
The dynamical exponent z, computed using the finite-size 
scaling of the gap, A oc L~ z , is shown in Fig. 2(b). We 
see that z = 2 c iean = 1 f° r D < D c re 0.3, whereas for 
stronger disorder, z becomes a monotonically increasing 
function of D. 

The data in Figs. 1 and 2 confirm our prediction that 
the ferromagnetic quantum phase transition is governed 
by the clean fixed point for weak correlated disorder. For 
stronger disorder, it is governed by a line of finite-disorder 
fixed points as shown in Fig. 1(c). This line of fixed points 
also emerges from an extension of the strong-disorder RG 
(SDRG) [9] to correlated disorder [19]. Moreover, for very 
strong disorder, the result z re D can be understood as 
arising from the singularity of the bare disorder distribu- 
tion via local fluctuations of weakly connected spin clus- 
ters. This line of fixed points belongs to a novel univer- 
sality class dominated by disorder effects that are pertur- 
batively irrelevant. In this class, some critical exponents 
(such as v) are universal and take their clean values, while 
others arc nonunivcrsal, such as z. 

Entanglement entropy — We now turn to the ground-state 
(GS) entanglement properties at criticality in the case of 
correlated disorder, Ji = hi. We first study the von Neu- 
mann entanglement entropy S(£) = — Tt{px In Pa) ; where 
Pa is the reduced density matrix of a subsystem A of size 
I. Because we wish to relate S(£) to local fluctuations of 
a globally conserved operator, we consider the spin-1/2 
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Fig. 2: (Color online) (a) Disorder averaged correlation length 
£ as a function of 8 (sizes up to 2 250 sites, averages over 20 
disorder realizations), (b) z vs. D at criticality 5 = (sizes 
up to 2 15 sites, typically 1 000 disorder realizations). In both 
cases, disorder is correlated, hi = e s Ji, and error bars are of 
size of the symbols. 
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h xx = aJ2u (s*s* +1 + sfs! +1 ) 



(3) 



from now on. It can be mapped onto the Hamiltonian 
(1) by setting t2i-i = hi and tn = Ji. The two systems 
share the same entanglement properties at criticality with 
25i sing (£) = Sxx{£) [20]. 

Figure 3 shows our numerical results for the average en- 
tanglement entropy S(£), computed via standard methods 
[13]. As expected, for weak correlated disorder, D < D c , 
S(£) ~ (c c ff/3)ln^, is universal and c G ff takes the clean 
value c e g = c = 1. Surprisingly, for D > D c , the en- 
tanglement entropy increases monotonically with D, i.e., 
S ~ jrc ffln£, with c c s > c, as shown in inset (a). In 
contrast, for uncorrelated disorder, c c s = cln2 < c [12]. 
The increase of S(£) with correlated disorder was first no- 
ticed in Ref. [21], but the different nature of the GS, its 
dependence on the disorder strength and its universality 
class were not considered. We also point out that a small 
deviation from perfect disorder correlations will drive the 
system back to the IRCP physics, as revealed by the SDRG 
and checked numerically [19]. 

Recently, it has been noted that S(£) can be related to 
the local fluctuations of certain thermodynamic quantities 
[22]. Here, we use this relation in order to gain further 
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Fig. 3: (Color online) Disorder averaged entanglement entropy 
S vs. subsystem size I for the critical GS of (3) with correlated 
disorder (hi = Ji) for various disorder strengths D (chains of 
1 024 sites averaged over 5 000 disorder realizations). Inset (a): 
effective central charge c c g vs. D. Inset (b): magnetization 
fluctuations amplitude n e g vs. Ceg • 

insight on the GS of (3). Hence, we compute the fluc- 
tuations T m if) of the magnetization 5j[ = J2%eA &i °f a 
subsystem A with length I: T m (i) = {{S% - {S Z A )) 2 ). In 
the clean case (U = t), the fluctuations were shown to be 
proportional to the von Neumann entropy, J- (I) = K In £, 
with k = 7r~ 2 , yielding c = tt 2 k [23]. For correlated dis- 
order (<2i-i = t2i), J 7 ^) = K f[hi£ where K e ff takes its 
clean value for D < D c and increases similarly to c c g for 
D > D c , as shown in inset (b) of Fig. 3. Interestingly, the 
ratio c c ff/n e ff — ir 2 remains identical to its clean value in 
the entire £>-rangc studied. Thus even though the line of 
critical points at D > D c is dominated by disorder, the 
relation between fluctuations and entanglement entropy 
strongly resembles the delocalized clean system. In con- 
trast, in the random-singlet state arising for uncorrelated 
disorder, it is easy to show that F(£) = K e g(hi£) with 
KeB = 1/12 [24]. 

Currently, to the best of our knowledge, there is 
only one general framework in which c c ff can be under- 
stood in terms of random-singlet phases of ID systems 
[25]. We remark that such framework does not apply to 
our correlatcd-disorder case as there is no corresponding 
random-singlet state. Thus, a fundamental understand- 
ing of c e ff here reported is still lacking and further funda- 
mental understanding of the entanglement entropy in ran- 
dom systems is desirable. Studing the correlated-disorder 
effects in the many random-singlet states considered in 
Ref. [25] , and comparing their entanglement entropy with 
their clean counterparts seems to be a helpful direction of 
research. 

Experimental realizations — Consider a quantum 5 = 1/2 
chain with nearest-neighbor exchange Js-s- Doping the 



chain with a small concentration of Si mp = 1/2 impuri- 
ties having Js ita -s Js-s would realize a slightly mod- 
ified SU(2) (Heisenberg) version of a correlated double- 
bond Hamiltonian (3), closely related to the Kondo prob- 
lem in spin chains [26, 27]. Based on the Matsubara- 
Matsuda representation of spin-f /2 degrees of freedom by 
hard-core bosons [28], one could also implement the cor- 
responding hard-core bosonic model (with now correlated 
random hoppings) using cold-atom systems. Finally, one- 
dimensional polymers such as the family of polyaniline 
[29] are modeled as random-dimer tight-binding chains 
[30] which can defy Anderson localization. This random- 
dimer chain can be mapped to the Hamiltonian (3). Our 
results thus provide an alternative view of the absence of 
localization in certain ID electronic systems. 
Discussion and conclusions — In summary, we have pre- 
sented a general mechanism by which local correlations 
between the random variables render a clean critical point 
stable against weak disorder even though it violates the 
inequality dv > 2. Although this appears to contradict 
the Harris criterion, we emphasize that it merely violates 
one of its preconditions, namely the spatial variation of 
the distance from criticality. Indeed, Harris [4] argued 
that a clean critical point is stable if the mean (local) dis- 
tance from criticality [r] av is larger than the width Ar 
of its distribution. For uncorrelated disorder, this yields 
Ar/[r] av ~ 2 < f for £ — > oo, recovering the (dv > 2)- 
form of the Harris criterion. However, for our correlated 
disorder, Ar = 0, and [r] av ^> Ar is satisfied regardless of 
the values of d and v. 

Our mechanism for suppressing the disorder effects by 
local correlations will also operate in d > 1. The ques- 
tion under what general conditions the random mass term 
can be removed completely will be relegated for future 
research. Interestingly, Yao et al. [31] recently reported 
apparent violations of the Harris criterion in several dis- 
ordered dimerized spin models. In one system (random 
dimer model), this has been attributed to the fact that 
the quantum critical point does not depend on the disorder 
strength, which strongly resembles our mechanism for the 
absence of random mass. A similar argument was given 
for the Mott-insulator to supcrfluid transition at the tip 
of the Mott lobes [32]. Another system of Ref. [31] (ran- 
dom plaqucttc model) shows a dependence of the critical 
coupling on the disorder strength, albeit a very weak one. 
We emphasize that in such a case, our mechanism would 
restrict the deviations from clean critical behavior to a 
narrow interval around the critical point that may well be 
unobservable. 
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